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ABSTRACT 



The Broclc version of the Split-Step Fast Fourier 
Transform (SSFFT) and the Jeager version of the Implicit 
Finite Difference (IFD) acoustic parabolic equation models 
are compared with a Lloyd mirror interference pattern in the 
range domain. The SSFFT displays the inability to place the 
transmission loss nulls at the correct ranges. It is also 
unable to utilize bottom loss information correctly. The 
IFD produced nulls at the correct ranges; however, it 
inserted an unacceptable amount of noise except when siaall 
(1 m) vertical grid steps were used and the pressure release 
bottom was placed at extended depths. In shallow water 
cases, the IFD is able to properly represent the pressure 
information. Each model is explored in the wavenumber 
domain by use of a "Wavenumber Technique" (WT) model with 
emphasis on source depth determination. The source depth 
may be determined by measuring the distance between the 
equally spaced nulls in the wavenumber representation. 
Neither acoustic model was able to provide accurate source 
depth information when the null spacings were compared to a 
known source-depth determination curve. Since the null 
spacings were not uniformly spaced, this was to be expected. 
Some specific proolem areas in the models were identified by 
the use of the WT. 
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I. IHTR0D3CTI0N 



Conventioaal aethods used to represent underwater sound 
transmission have depicted tra nsmissi on loss (d3) as a func- 
tion of range in either tabular or grapnical form. Since 
the interest was in determining quantities such as range of 
bottom bounce, convergence zones, and probability of 
detection, these methods proved very satisfactory. However 
when it became necessary tc compare various underwater 
acoustic propagation loss models with each other or actual 
in-s itu experiments, they provided the analyst with less 
than adequate insight into the causes for the observed 
errors in the various models outputs. Furthermore, interest 
has been generated in analyzing the received signal for 
information from the sound source which cannot easily be 
induced from transmission loss as a function of range. 

One method, the Ravenumber Technique (WT) , described by 
F. R. DiNapoli of N3SC £Ref. 1] and applied by Richard Lauer 
of NCRDA £Eef. 2h provides more information for model 
comparisons and may have the capability of determining sound 
source depth and range from the receiver in certain cases. 
DiNapoli used an analysis of acoustic propagation in the 
wavenumber domain as an internedxate step in getting the 
transmission loss curve in the Fast Field Program (FFP). 
Lauer studied the wavenumber domain information from the FFP 
and concluded that source localization mignt ce possible. 
B. B, Stamey [Ref. 3] investigated the potential for deter- 
mining the source depth by using the Brock Split-Step Fast 
Fourier Transform (SSFFT) [Ref. 4] parabolic equation model 
with the "Ravenumber Technique" (WT) model. His preliminary 
investigation involved the use of isospeed and ASTHEX sound 
speed profiles with a fully ansorbing bottom, varied 
source/receiver combinations, and multiple frequencies. 



The present investigation uses two acoustic jarabolic 
equation models for comparison with theoretical Lloyd mirror 
depictions, and the effects of observed inconsis tencies in 
predicted transmission loss curves on the calculated results 
from application of the WT model. An analysis of the sensi- 
tivity in determining source depth based on model output is 
attempted to ascertain possible use in naval operations. 

To facilitate a straight-forward analysis of the results 
and evaluation of the HT technique for the cases studied, 
only a selected number of the various computer runs are 
included. Only those curves which best illustrate the 
effects revealed oy this research are presented. 



II. W A7EH0H BER TECHNIQUE (WT) 



A. BACKGEODND (LLOYD MIEBOfi) 

The WT can be elucidated with the help of the classical 
Lloyd mirror effect [Eef. 5] which describes the interaction 
between direct path and surface reflected sound signals from 
the same continuous wave (CW) source. Figure 2. 1 illus- 
trates the geometry involved in the Lloyd mirror effect. 




Figure 2.1 Lloyd Mirror Geometry 



The ccmplex pressure (P) can be expressed by 



P= A 



1 -jkr; 1 -jkrc 



fi 



(eqn 2.1) 






when r: and are related by 



= VR"^(d-h)"’ r, = ^JWT{d7h? 



(egn 2.2) 



with the assumption r^ the pressure amplitude (p) can be 

approximated by 



p _ 2fhd m =1 2 3 



(egn 2.3) 



The interaction between the direct path and surface 
reflected waves produces constructive and destructive inter- 
ference which is manifested as peaks and nulls for the in 
phase and out of phase conditions respectively. Figure 2.2 
was produced by using eguation 2.3. The null locations can 
be obtained by the relationship. 



P = 




(egn 2.4) 



Shen eguation 2.3 is expressed in the horizontal wavenumber 
domain (k^)/ the details of the resultant functional depen- 
dency on horizontal wavenumber yield useful information. 
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Figure 2.2 Theoretical Lloyd Mirror Transmission Loss. 



B, PHYSICAL PROCESS 



The WT is a process by 
(with range dependent ampl 
into the spectral density 
wavenumber. The use of th 
ment would require a guadr 
signal to attain the com 
demodulation process can be 
in figure 2.3 [Ref. 6]. 



which 


. the com 


plex pr 


essur 


e w 


ave 


itude 


and phase 


) is 


tr ans 


for 


med 


as a 


function 


of the 


hori 


zon 


tal 


e HT 


in the op 


eration 


al en 


vir 


on- 


ature 


demodulat 


ion of 


the 


sou 


rce 


plex 


pressure. 


The 


quadrat 


are 


St be 


defined b 


y the illust 


rat 


ion 



2cos(0J.t) 






Figure 2 


,3 Qu 


Rhen 


acoustic 


models 


the 


requirement to 


eliminated 


since 



Quadrature Demodulation of a Signal, 



perform a quadrature demodulation is 
the complex pressure is directly 



16 



accessible. In actuality, acoustic parabolic equation 
models provide the complex envelope function (U) which must 
be multiplied by the HanXel function in order to obtain the 
complex pressure. The acoustic wave is written in the form 

P = UHo’(k 2.5) 

where the reference wavenumber is defined by 

ko = -^ (egn 2.6) 

V O 

and 



U = U(r,z) 



(eqn 2.7) 



is the solution to the appc 
complex pressure is corre 
Fourier transformed to atta 
displaying the results grap 
number on the x-axis and t 
the y-axis is constructed a 
measured. Figure 2.4 [from 
for a source and receiver i 
Since the Lloyd mirror 
given by [Ref. 2] 



opriate parabolic equation. Ihe 
cted for volume attenuation and 
in tne spectral density. A plot 
hicaliy with the horizontal wave- 
he normalized spectral density on 
nd the spacings between the nulls 
Ref. 2 ] is an example of the WT 
n the same type of water mass, 
field in tae wavenumber domain is 



F(k) = 



s'\r.U3Zt) 

P 



ej</32r) 



(eqn 2.8) 
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an alternative formulation [Ref. 2], which produces evenly 
spaced nulls, is obtained by converting the horizontal wave- 
number to beta 




(eqn 2.9) 



By using 

distance, 

ascertain 



beta instead of the 
delta beta, between 
the source depth by 



horizontal wavenumber, the 
any two nulls can be used to 



A/? 



- JL. 

■ z. 



(egn 2,10) 



as illustrated in figure 2. 

The application for im 
isospeed cases, to describ 
pressure release boundary 
Therefore the parabolic egu 
Lloyd's mirror effect and 
issue, cf. Chapter 4 of 3r 

Fhen the WI is impleme 
ensure that the complex pre 
of points to describe the 
sure field. This will pre 
in the wavenumber domain, 
of the pressure signal mus 
begins and ends with zero 
necessary because tne pr 
repetitive oscillation for 
sure array is zero-filled 
generate an array which ha 
the Fast Fourier Transform 
increment is generated with 



5. 

plying the method of images, in 
e the erfects of a perfectly fiat 
can be justified by inspection, 
ation can be used to produce the 
surface decoupling is not an 
ekhovskikh [Ref. 7]. 
nted, care must be exercised to 
ssure contains an adequate number 
shortest periodicity in the pres- 
clude the possibility of aliasing 
The real and imaginary elements 
t be modified so that the signal 
values. This modification is 
ensure signal must represent a 
the transform. The complex pres- 
beyond the data in order to 
s the length of a power of 2 for 
(FFT) . The Horizontal wavenumber 
in tne computer code by 



18 



(eqn 2. 11) 



krm 



= ko-^ 



27 T / m 
Ar NPT 
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Figure 2.5 Source Depth Determination Curve 



III. A C0D5 TIC MODEL AN ALYS IS 



A. PARABOLIC EQUATION 

The elliptical wave eguatioc can be approximated by a 
parabolic equation (PE) when it is assumed that the envelope 
function varies slowly with range. A detailed mathematical 
description of the PE from an acoustic point of view is 
presented by DeSanto [Ref. 8 and 9] and a simpler but less 
general description can be found in Coppens [Ref. 10]. The 
parabolic equation has the form 



where it is assumed that the pressure has the form of equa- 
tion 3.2, 



S(r) represents the primary radial dependence of the field 
in terms of an outward propagating cylindrical wave of the 
form [Ref. 4] 



If it is assumed that the range of interest is many wave- 
lengths from the source; then the asymptotic form of the 
Hankel function (equation 3.4) can be used [Ref. 11]» 




d Z 



P = U(r,z)S(r) 



(eqn 3.2) 



S( r ) = n'P(kor) 



(eqn 3.3) 
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Two PE models# the Brock version of the Split-Step Fast 
Fourier Transform (SSFFT) and the deader version of the 
Implicit Finite Difference (IFD) were investigated for 
compatibility with the HT in determination of source depth. 

B. SSFFT 

1. Transm i ssion loss C omp arison wit h Lloyd M irror 



needed a model to provide the pressure as a function of 
range. Initially the SSFFT model [Ref. 4] was selected . 
This model generates successive values for U as a function 
of range with the help of the algorithm 



This model was a natural choice since it was the model used 
in the preliminary study dealing with the determination of 
the source depth by the WT [Ref. 3]. However# it soon 
became apparent that this model’s inherent weaknesses would 
require the consideration of another model if ocean bottom 
interactions were to be studied. These weaknesses# and 

their impact upon the MT# will be discussed shortly. 

The SSFFT is a range dependent acoustic wave model 
which# for this analysis# will be operated in a range inde- 
pendent manner. In other words# only one sound speed 

profile will be used, the water mass will be assumed homoge- 
neous# and the bottom flat. To facilitate the study of the 
SSFFT model# the source code of the model was modified to 
generate the variables required for follow-on programs to 



To analyze sound waves in the wavenumber domain# we 







(eqn 3.5) 
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transforn the pressure information and then display the 
results graphically. A copy of the source code listing for 
each program is provided in appendices A, 3, and C. These 
programs were then linhed by Job Control Language (JCI) so 
that the WT graphics were an automatically executed. 

A test run of tne SSFFI was made for 100 hz, 500 m 
source and receiver depths, 0.019 km range step, fully 
absorbing bottom, 2000 m water depth, and an isospeed 
profile at 1500 m/sec. Tnese values were chosen so that a 
comparison could be made with the tnaoretical Lloyd mirror 
transmission loss (figure 2.2). Ine output was displayed as 
a transmission loss curve (figure 3.1). From equation 2.4, 
the location of null number 1 (R1) should be 33.33 km and 
null number 2 (E2) should be at 16.67 km. The SSrFT placed 
R1 at approximately 38.3 km and R2 at 17. 1 km. Since the 
model kept the frequency, source depth, and receiver depth 
constant, then this equates to a reference sound speed of 
1305 and 1462 c/sec respectively. It is apparent that as 
the range decreases the error is reduced which inaicatss the 
possibility of a problem in the range step section of the 
model which leads tc the introduction of a systematic and 
increasing deviation of the calculated sound field from that 
predicted by the classical Lloyd’s mirror interference 
pattern. The "washing out" of the interference pattern for 
ranges less than about 2 km is a result of the particular 
approximation of a point source used to initiate the 
program, which is tc be expected. Ratner regular fluctua- 
tions develop beyond about 30 km; from the geometry of the 
case, it is plausible that this arises from interference 
with bottom reflected signals. 
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Figure 3.1 Test of SSFFT Model for Null Spacings. 



Another test was made using the same inputs as 
that tne source and receiver were moved to a dep 
(figure 3.2). Since the half beaawidth is 23 
the bottom is fully absorbing, the snortest rang 
surface reflection cculd occur is approximat 
beyond this range, surface reflections should 
accordance with equation 2.4. The nulls beyon 
occur at the correct ranges; however, there are 
tions of the overall transmission loss signal, 
modulations are another indication that 



fully absorbing 
program. 
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Figure 3.2 Test of SSFFT Model for Bottom Interaction 
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• iil Ana l V sis 

Ihe wr cutout from 
units, wiiich is a problem 
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expressed by equation 2.6, 
hz with a reference sound 
Q/sec) is 0. 1277 1/ft. A 
wavenumber for a range in 
points used to represent 
maximum beta value shoul 
reasonably well with the ma 
3.4, ana 3.5. Contrary to 
spaced and the expected d 
(equation 2.10) is net rep 
ings. Stamey £2ef. 3], 
spacing of beta nulls. Th 
not represented by any of t 
tne improper null placement 
of equal spacing between 
something else, possibly 
represented as a continuous 
within the source code, 
difference, it must be reso 
with tne KT to determine th 



the S3FFT model is in English 
when comparing different models, 
tion the problem was minimal. As 
the reference wavenumber for 100 
speed of 4921.26 ft/sec (1500 
fter correction of the reference 
crement of 62.3 ft and the 2048 
the wavenumber spectrum, the 
d be 0.1248 1/ft which agrees 
ximum beta values in figures 3.3, 
theory tne nulls are not equally 
elta beta value of 0.00192 1/ft 
resented by any of the null spac- 
observad cne same non-uniform 
e reason why the source depth is 
he null spacings may be traced to 
in the range domain. The lack 
the nulls must be attributed to 
the pressure wave was improperly 
signal or numerical errors exist 
Whatever the reason for this 
Ived before the S3FFT can be used 
e source depth. 
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Figure 3.4 SSFFT WT Plot with Receiver at 1000 ueters 
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Figure 3.6 IFD Environment. 
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Figure 3.7 IFD Transmission Loss. 
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Figure 3.8 IFD Detailed Transmission Loss 
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Figure 3.9 IFD Noise Free Transmission Loss 
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Figure 3.10 IFD: Shallow Water Transmission Loss 
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HT plots after being Fourier transforned. Continued anal- 
ysis revealed that the noise could be elininated if the 
thicknesses of the artificial attenuation layer and the 
bottom pressure release surface are chosen so that they are 
at least twice the water depth. The choice of a small 
vertical grid step will reduce the thickness required for 
the pressure release surface, but this version of the IFD 
program is limited to 5000 vertical grid points due to soft- 
ware restrictions. Figure 3.9 is a transmission loss curve 
obtained from a 500 m source/receiver depth, a 2000 m water 
depth, a 4000 m upper level artificial attenuation layer, 
and an 8000 m bottom pressure release surface. Comparison 
of figures 3.1 and 3.9 with the classical Lloyd's mirror 
transmission loss (figure 2.2) snows that tne shape of the 
curve beyond 30 km is greatly improved with the IFD. The 
number of vertical grid steps was 5000. When the IFD is 
used in a shallow water situation and at 25 hz, the noise 
does not appear to be present. Figure 3.10 is a plot of the 
transmission loss for a signal at 25 hz, witn the source/ 
receiver depth at 50 a, the bottom sloping up from 350 to 50 
m, and an artificial attenuation layer and bottom pressure 
release surface located at 750 and 1000 m, respectively. 
The noise whizh was so evident in figure 5.7 is absent from 
the shallow water case. Producing tne noise free curves 
requires excessive computer time except for shallow water 
cases. 
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2 . WT Ana ly sis 

To reduce the amount of the pressure information and 
yet preserve the details, the IFD tests were conducted 
between 10 and 40 Jem in range. The reference wavenumber for 
tne test input with the IFD was 0.4189 1/m which, after 
correcting for a range increment of 7.5 m with 4096 points, 
produced a beta maximum of 0.4188 1/m. Figures 3.11, 3.12, 
and 3.13 indicate agreement with the beta maximum and a 
delta beta of 0.00628 1/m can be found among the null spac- 
ings. When the vertical grid step is reduced in size, the 
output signal contains more information of finer detail. 
The finer detailed information will, if noise is present, 
produce more clutter on the WT depiction (figure 3.13). In 
numerous WT depictions a peak was observed at the minimum 
and maximum wavenumbers. Stamey believed that the right and 
left intensity maxima corresponded to beam elevation angles 
of 0 and 30 degrees, respectively. He further stated that 
the left maxima also may be related to algorithmic inaccura- 
cies [Ref. 3]. The "U-shaped phenomenon" observed during 
Stamey' s investigation of the S5FFT, is present with the IF3 
and the nulls are not equally spaced. 
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Figure 3.11 IFD: 500 meter Source, 15 meter Vertical Grid Step 
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Figure 3.12 IFD: 500 meter Source, 3.75 meter Vertical Grid Step. 
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Figure 3.13 IFD: 500 meter Source, 0.8 meter Vertical Grid Step 



ffne reason for the "O-'Shaped phenomenon” appears to 



related to noise in the envelope function field which is 
further amplified when the Hankei function is inserted to 
obtain the complex pressure. Using a low pass filter and 
increasing the vertical grid step 
problem. In figure 3.8, the null 
km when a single null should exist 
tion or noise produces the mos 
masking the desired density spectrum, 
filter would probably solve this 
but eventually the IPD will regui 
nate the noise. '/?hen the shalic 
(figure 3.14), the "U-shaped pher 
null locations are more distinct, 
water case, the bottom interact! 
interference but the null spacin^s 
value for the source depth (0. 0623 1/m). 
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17. C 0MCL 3SI0NS 



Two paratolic equation models, the SSFFT and IFD, were 
used to predict sound fields for comparison with the Lloyd’s 
mirror interference pattern in the range domain and then in 
the wavenumber domain. In the range domain, the SSFFT 
improperly handled the location of pressure nulls and 
displayed bottom interference. While these weaknesses in 
the SSFFT were absent in the IFD, the inadvertent and dele- 
terious insertion of noise was extensive. In tne wavenumber 
domain, the noise in the IFD was prohibitive in all but the 
deep sediment and shallow water cases. The IFD produced 
results which were comparable with the SSFFI when the noise 
was not a factor. The range step and tne bottom attenuation 
algorithms of the SSFFT should be reviewed in order to 
correct those problems described above. In the IFD, the 
noise problem will require further investigation. If the 
problem in tne IFD is not corrected, this model will be 
severely limited. With the guidance provided by Lauer 
[Ref. 2], the wavenumber domain can provide the detailed 
information needed to quickly correct inconsist encies 
observed between acoustic models. Therefore, the FT should 
be investigated further to ensure that any software weak- 
nesses are eliminated and a thorougn understanding of the 
process is assured. 

Experiments in the operational environment should be 
planned to test the NT with actual sound sources. The ship- 
board experiment could be easily accomplisned since the NT 
can be implemented on micro-computers if the pressure infor- 
mation is provided from an outside source. The NT, if it is 
going to be considered as a method to determine source 
depth, appears to have three major snortcomin gs: 
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1. A clean CV signal with pnase infornation is essen- 
tial, which will reguire quadrature demodulation and 
extensive filtering. 

2. Prediction models are required which will provide 
equally spaced beta nulls so that software can be 
developed to automate the process. 

3. High resolution processing equipment is necessary 
because the change in beta spacing for shallow sound 
sources is small (figure 2.5). 
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APPENDIX A: SSFFT SOURCE CODE (AT NPS) 
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IPRNT = 0 DO NOT PRINT TRANSMISSION LOSS TABLE 

.NE.O PRINT OUTPUT TRANSMISSION LOSS TABLE 
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HETUHN CODE - FLAG = 0.0 SUCCESSFUL ADVANCE 
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DEFINE THE iSVEEAGl ACOUSTIC WAVELENGTH AND WAVE NUMBER 
AND THE MESH INCREMENT IN TRANSFORM SPACE. 
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GET VOLUME ATTENUATION FACTOR 
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INPUT - F5CTOR = DR ♦ K / 2 
FN = N**2 - 1 

ATTEN = VOLUME ATTENUATION 
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INTRODUCE ATTENUATING MUD SEGMENT TO PROFILE 



FNO=FN (HW) 

IF (NBOIM.L‘1.0) FN 1 = FNMUD(FN0, 0) 
Fl=j3. 1415S26535/HORBAN) 
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OUTPUT - TICSS TRANSMISSION LOSS AT DEPTH 
FUNCTION TICSS (RR,Z) 
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INTERPOLATE N**2 - 1 ON THE FIELD HESH. 



ZW=DZ*FL0A1 (NH) 

IF (HC.EQ. 1 ) ZW=ZW*SQET (CO/SPEED (ZW) ) 
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END 

SUBROUTINE R8SYN ( INI , DO , B 1 , 32 , D3 , D4 , B5 , B6 , E7) 
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THE FUNCTION ALPHA IS CONSTPUCTEC AT MESH POINTS GIVEN 
BY ZHAT().. ALPHA IS BEPRESENTEC BY A PIECEWISE 
CONTINUOUS FUNCTION IN DEPTH, WITH THE COEFFICIENTS 
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CVP = A'lIN 1 (CVP,C (J) ) 
CCNT INUE 
WRITE ( 2 ) CVP 
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APPENDIX D: IFD SOURCE CODE (AT NPS) 
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BOTTOM SLOPES DOWN SLOWLY, BOTTOM MODIFIED 
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IS SECTION PERTAINS TO POINTS IN KATER COLUMN 
10 I=1,NWHAX 

*** CALCULATE REAL INDEX Oi- REFRACTION IN EATER 
XN = CO/CWATER(I) 

*** CALCULATE ATTENUATION AS PER COMMENTS IN SU2R0UTIN 
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INITIAIIZE RANGE VARIABLE AT WHICH SOLUTICN IS TO BE RECORDED 
XHF = EAU-WDR 



♦** COMPUTE RECEIVER DEPTH TO NEAREST DZ 
IF { 2R.LT.DZ ) ZH = DZ 
IWZ = 2E/DZ + 0.5 
ZR = IV2*DZ 
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UPDATE X MAIHIX ELEMENTS ON INTERIACE 
C(IfACE2,1) = XLI 

C’IBACE2,2) = -AJIFACE2) * ZZ5 + ZZIC 
C IfACE2,3) = XRI 
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(1) THIS SUBROUTINE IS EFFECTIVELY THE COHTINUAIION OF 
SUBROUTINE PRTHTl. 

(2) THE FILE CREATED CORRESPONDS TO UNIT FILE NUdDEK; 
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